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ABSTRACT 

We analyze the redshift evolution of the Rotation Measure (RM) in iTavlor et al.l (|2009| ) dataset, 
which is based on NVSS radio data at 21 cm, and compare with results from our previous 
work (jKronberg et a l. 2008; Ber net et al.ll2008l 120101 ) . based on RMs determined at lower wavelengths, 
e.g. 6 cm. We find that, in spite of the same an alysis, Taylor's dataset produces neither an increase 
of the RM dispersion with redshift as found in iKronberg et ail (120 08) . nor the correlation of RM 
strength with Mgll absorption lines found in iBernet et al.l (|2008f h We develop a simple model to 
understand the discrepancy. The model assumes that the Faraday Rotators, namely the QSO's host 
galaxy and the intervening Mgll host galaxies along the line of sight, contain partially inhomogeneous 
RM screens. We find that this leads to an increasing depolarization towards longer wavelengths and 
to wavelength dependent RM values. In particular, due to cosmological redshift, observations at fixed 
wavelength of sources at different redshift are affected differently by depolarization and are sensitive 
to different Faraday active components. For example, at 21 cm the polarized signal is averaged out 
by inhomogeneous Faraday screens and the measured RM mostly reflects the Milky Way contribu- 
tions for low redshift QSOs, while polarization is relatively unaffected for high redshift QSOs. Similar 
effects are produced by intervening galaxies acting as inhomogeneous screens. Finally, we assess the 
performance of Rotation Measure synthesis on our synthetic models and conclude that the study of 
magnetic fields in galaxies as a function of cosmic time will benefit considerably from the applica- 
tion of such a technique, provided enough instrumental bandwidth. For this purpose, high frequency 
channels appear preferable but not strictly necessary. 

Subject headings: galaxies: high-redshift — galaxies: magnetic fields — quasars: absorption lines — 
galaxies: evolution 



1. INTRODUCTION 

Faraday Rotation Measures (RM) is one of the very 
few methods to probe extragalactic magnetic fields. The 
RM is given by the change in observed polarization angle, 
Axoj over a change in the observed wavelength square, 
AAq. For a polarized radio source at cosmological red- 
shift z s it is defined as 

o 

„, f , x Axo ai 1fl5 [ n e (z)B {{ [z) dl 

where the RM is in units of rad m~ 2 , the free elec- 
tron number density, n e , is in cm -3 , the magnetic field 
component along the line of sight, B\\ , is in Gauss, and 
the comoving path increment per unit redshift, dl/dz, is 
in pc. Eq. [T] assumes a uniform RM screen across the 
source and a spatial separation of the linearly polarized 
sour ce and the Faraday rot ating plasma. 

In IKronberg etaLl (f2008h (K08) we used a sample of 
268 RM values of extragalactic radio sources to assess 
the redshift evolution in the RM dispersion. We found 
an increase in the RM dispersion with redshift, which 
became statistically significant above z ~ 1. We pos- 
tulated that the increase in the RM dispersion is pro- 
duced by magnetic fields in intervening galaxies. To test 
this hypothesis spectra of 71 QSOs at UVES/VLT were 
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taken. In IBernet et al.l (|2008t ) (B08) we showed that in- 
deed sightlines with intervening strong Mgll absorption 
systems have significantly higher RM values than those 
without. The findings in that work implies that ~ 10 /iG 
magnetic fi elds exist in galaxie s out to z ~ 1.3. 

Recently ITavlor et all (pool (TSS09) determined RM 
values of 37'543 sources, based on polarization observa- 
tion of the NRAO VLA Sky Survey (NVSS). After the 
release of the RM catalogue we compared our RM values 
with those of TSS09 and found major differences between 
the two datasets. As we show in section 12.21 using the 
RM values of TSS09 the results of K08 and B08 can not 
be reproduced. In this work we present an analysis of 
the differences found between the two datasets. We then 
develop a toy model to show that such differences can 
be accounted for by inhomogeneities in the RM screens 
of the QSOs host galaxies and intervening galaxies. In 
view of this model, due to the strong depolarization ef- 
fects, RM data based on low frequencies observations 
and interpreted according to Eq. (JTJ) are inadequate to 
probe inhomogeneous RM screens produced by interven- 
ing galaxies. 

Recently there has been an increased attention to the 
effects of inhomogeneous RM screens on the observed 
degree of polari zation and pola r izatio n angle as function 
of wa velength. iRossetti et al.1 (|2008| k iMantovani et al.l 
( 2009) used multiwavelength radio observations to model 
depolarization of their sources. They showed that the 
depolarization as a function of wavelength can be better 
described including a covering factor for the inhomoge- 
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neous screens. iFarnsworth et al.l ()201lD used Westerbork 
Synthesis Radio Telescope (WSRT) observation at ~ lm 
combined with NVSS observation at 21 cm to test differ- 
ent depolarization models and RM determination meth- 
ods. T hey consider ed the traditional Faraday dispersion 
screen ([Burn! [T966) . a two component model for Fara- 
day Rotation which produces oscillation in the degree of 
polarization, and the RM Synthesis method. The com- 
parison of the different RM values obtained from these 
methods revealed that if Faraday structure is present, 
the different methods may lead to different RMs. They 
further stress the importance of considering both the po- 
larization angle and amplitude for the correct determi- 
nation of the RMs. 

While several effects may be at work, as briefly dis- 
cussed in section [231 in this paper we keep the level of 
sophistication of the model to a minimum and focus on 
the role of Faraday screen particularly on simple ideas 
associated with inhomogeneity which can account for im- 
portant observational effects and test these ideas for con- 
sistency with available data. The rest of this paper is 
organized as follows. In section [2] the RM datasets of 
B08 and K08 arc compared and the differences are pre- 
sented. Previous work on depolarization of extragalactic 
sources is shortly presented in section [3J In section 2] a 
depolarization toy model is presented which can account 
for the differences in the RM values of the two datasets. 
A discussion of the nature of the partial inhomogeneous 
RM screens and computed Faraday spectra as compari- 
son for RM surveys is presented in section [5j In section 
[5] we give a summary of our findings. 

2. ROTATION MEASURE DATA 
2.1. Datasets 

The RM values used in K08, B08 and iBernet et all 
(|2010l) were collected by P. Kronberg and collabora- 
tors over the past decades using polarization observa- 
tions at various telescopes, including VLA and Effels- 
berg. The sample consists of 901 sources with de- 
termined RM and redshifts, partly taken by Kron- 
berg and collaborators, partly take n from the litera- 
ture ( Simard- Normandin et al.l 119811) . Until the work 
of TSS09 it has been the only existing large RM cata- 
logue. Generally at least three wavelengths were used 
for the RM determination (Phil Kronberg, private com- 
munication) using polarization data between 37 to 0.9 
GHz, with the bulk of the data betw een ~ 10.5 to 1.6 
GHz (jSimard-Normandin et "aTlll98lD . 

In the following sections we will assume that the deter- 
minations of the RM values of K08 were typically done at 
6 cm and designate the RM values as RMq. Of course due 
to the heterogeneity of the sample the wavelength range 
and typical value over which the RMs are determined 
might vary. While this is undesirable, its main effect 
should be the introduction of noise in the relations pre- 
dicted by our model between measured quantities. There 
might be additional issues associated with the heteroge- 
nous character of the polarization data used for the RM 
determinations, which, however, are not addressed here. 

The sample of K08 consists of 268 sources at Galactic 
latitudes |&| > 45° (exact definition given in K08). In 
B08 we obtained high resolution spectra of 71 relatively 
bright QSOs and relaxed the selection to |6| > 30°. Here 



we use the subset t hat w as employed in the study of B08 
and IBernet et all (|2010fl (the RM data are still propri- 
etary and will be published elsewhere by P. Kronberg), 
except when we look for differences in the redshift evo- 
luti on and use al l sourc es at |6| > 30° which were also in 
the lTavlor et al.1 ()2009f ) catalogue. 

The RM values of TSS09 are based on th e NRA O VLA 
Sky Survey (NVSS) from ICondon et all (fl99l which 
covers the sky at declinations |<5| > 40° in Stokes I, Q 
and U. The survey imaged the sky at 21 cm with a reso- 
lution of ~ 45 arcsecs and produced a catalogue of 2 x 10 6 
discrete sources. TSS09 choose a subsample of this cat- 
alogue and derived RM values of these sources based on 
determination of the polarization angles at 1364.9 MHz 
and 1435.1 MHz. The subsample was selected by request- 
ing a source intensity I > 5 mJy and a 8a detection in 
polarized intensity. To ensure that the polarized inten- 
sity is not dominated by instrumental effects they only 
considered sources with a fractional polarization greater 
than 0.5 %, which yielded a RM catalogue of 37'543 ob- 
jects. In the following sections we will designate the RM 
values of TSS09 as RM 21 . 

2.2. Discrepancies 

After the release of the TSS09's RM catalogue we 
checked if it contained any of the sources in B08, and 
found this to be the case for 54 out of 71 of them. Some- 
what surprisingly we found striking differences in the RM 
values which are briefly presented here: 

i) A comparison of the RM values from K08 and TSS09 
reveals that the biggest differences in the RM datasets 
correspond to large RM values of K08 and sources having 
a low degree of polarization p2i< Here p2i is the degree of 
polarization at 21 cm from TSS09 (average degree of po- 
larization at the two wavelengths 20.89 cm & 21.96 cm). 
The situation is illustrated in Figure [TJ Sources with 
P21 > 3.2% (black cross, 25 QSOs) have quite similar RM 
values in both datasets and mostly crowd around the di- 
agonal line. On the other hand sources with P21 < 3.2% 
show a big scatter around the diagonal line. Further- 
more there are sources which have large RM values ( 
|RM| > 50 rad uT 2 ) in the dataset of K08 but very 
low |RM| values in the sample of TSS09. The value of 
P21 = 3.2% maximizes these effects in Figured] 

ii) Using a sample of 371 RM values from TSS09 at 
latitudes \b\ > 30° and for which the QSO's redshifts, z, 
were available, we found no increase in the RM disper- 
sion with z, contrary to the main finding of K08. The 
situation is illustrated in Figure [5J where the cumula- 
tive RM distributions at 6 cm (left panel) and 21 cm 
(right panel) are shown, split at the median fractional 
polarization of the total sample, P21 = 2.5% (blue and 
black lines), and according to whether the source red- 
shift is below or above z = 1.0 (solid and dash line, re- 
spectively). The cumulative distribution of |RM 2 i| shows 
no evolution with z in both low and high p2\ samples, 
whereas there is a significant broadening of the low p2\ 
(0.5% < p 2 \ < 2.5%) |RM 6 | distribution with z. A 
Kolmogorov-Smirnov test (KS test) reveals that the RM 
distributions above and below z = 1.0 are different at 
a significance level of 99.75%. We also show (red dis- 
tributions) sou rces with determine d RMg values which 
are not in the I Taylor et al.l (|2009f) catalogue, but have 
declination \S\ > 40° and are thus in the NVSS. Most 
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likely these sources are below the inclusion limits of the 
iTavlor et~a l. (2009) catalogue, which are a 8a detection 
in polarized intensity and P21 > 0.5%. A redshift evo- 
lution can be seen also for these sources with a similar 
significance level of 99.64%. Most likely the broadening 
of the RM distribution towards lower polarization P21 is 
due to Galactic depolarization and possibly partly due 
to an increased error in RM, <5rm> towards lower po- 
larized flux levels. This is supported by the fact that 
Sum oc 5p/P, where P is the pola rized flux and dp its 
error (|Brentjens fc de Bruvnll2005H . and the high corre- 
lation between P21 and polarized flux, as attested by a 
Kendall's r test, which yields r = 0.40 and a chance 
probability of only 4.4 x 10~ 30 using all 371 RM values. 
However the broadening with z seen in the |RMg| distri- 
butions is hard to explain by an observational bias. The 
interpretation of a P21 dependent redshift evolution of 
the |RMg| distribution is however not straightforward ei- 
ther. In section 2] we show that it is crucial to distinguish 
between sight lines with/ without intervening galaxies. 

iii) The result of B08 that the RM distribution for 
QSOs with strong Mgll systems along their sight lines 
is broader (at ~ 95% significance level) than for QSOs 
free of absorbers, completely disappears using the RM 
values of TSS09. Figure |3] shows the cumulative distri- 
bution (dash histograms) of |RM 2 i| values from TSS09 
for sight lines with (iVMgii > 0, black)) and without Mgll 
(A^Mgii = 0, red) absorption systems. Clearly there is no 
difference between the two RM distributions. For com- 
parison the solid histograms show the distributions of the 
|RM 6 | values of B08. 

2.3. Caveats 

In the following sections we present a simple model 
to describe the impact of inhomogeneous, redshift de- 
pendent Faraday screens, on the polarization and RM of 
distant QSOs, as a function of wavelength. 

However, additional effects may be at work when com- 
paring quantities measured at different frequencies. In 
particular besides potential large uncertainties measured 
with both datasets, the low and high frequency emissions 
may originate from separate components in the source 
(e.g. the radio lobes and the compact core in a radio 
galaxy), inducing all sorts of frequency dependent effects, 
e.g., location, angular size, spectral indices and fractional 
polarization. For example, in some radio sources the frac- 
tional polarizatio n appear to increase and then decrease 
with wavelength (|Conwav et al.lfl~974f) . 

Also, while our model could be validated (or ruled out) 
by comparison of its predictions with observational data, 
this would require a high quality dataset, where the ef- 
fects discussed above (and possibly others) are kept un- 
der strict control. Unfortunately such a dataset is not 
available to us at the moment, which forces us to refrain 
from such comparison. The latter, however, might be 
possible in the near future, with the delivery of new ra- 
dio data from polarization surveys, either planned or just 
under way (see Sec. I5.2j) . 

3. DEPOLARIZATION BY INHOMOGENEOUS 
FARADAY SCREENS 

3.1. Previous work 
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Fig. 1.— Comparison of the RM values of K08 and TSS09, RM 6 
and RM21 respectively, splitted according to whether the degree 
of polarization is above (black crosses) or below (red filled circles) 
P21 = 3.2%. The two datasets differ especially for low P21 values. 
Errors on RMg are available for 24 of 54 RM measurments. 



In this section we develop a simple model to under- 
stand the effects arising from inhomogeneous Faraday 
screens and show that these can produce important de- 
polarization effects that reproduce the differences of the 
RM datasets of K08 and TSS09. 

The observed polarization p can be written as a com- 
plex number as: 

P = pe 2lx , (2) 

where p = P/ 1 is the fractional degree of polarization 
given by the ratio of the polarized and total intensity, 
P and I, respectively and x is the polarization angle de- 
termined by the other Stokes parameters, U and Q, as 
X = l/2arctan((7/Q). 

The linear dependence between the polarization an- 
gle, x an d A 2 stated in Eq. Q] is only valid for the 
case of a uniform foreground screen. Both the pres- 
ence of unresolved inhomogeneities in the Faraday screen 
and/or sources of polarized radiation embedded within 
the Faraday active region produce increasi ng depolar 



ization of the source at longer w avelength (|Bur: 



lepolar- 
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lTribbld[T99ltlSokoloff et al.lfl998l) . This leads to nonlin- 
ear dependencies of the polarization angles and A 2 . In 
such cases detailed modeling is necessary for the correct 
interpretation of the RMs and for their use in measure- 
ments of magnetic fields. 

There are several mechanisms that can reduce the frac- 
tional degree of polarization either in the radio source 
itself or in its foreground: 



i) iBuTnl (fl966T) : ISokoloff et all (Il998l) showed that 
random fluctuations in the magnetic field within 
the source lead to a wavelength independent re- 
duction of the degree of polarization. The latter is 
given by the ratio of the regular-to-total magnetic 
field energy, p = p t (l)B 2 / B 2 , where p t (-y) = f^pf 
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Fig. 2. — Comparison of the cumulative distributions of RM values measured typically at 6cm (left panel) and at 21 cm (right panel). The 
samples are split according to the polarized fraction P21 and the redshifts of the sources. Significant redshift evolution is seen in the |RM6 
values but not for the |RM2i| values. For sources with 0.5% < P21 < 2.5% the RM distributions for sources above z=l is broad er than for 
sources below z=l, at a significance level of 99.75%. Red histograms are for NVSS sources not included in Taylor et al. (2009) catalogue. 
This is because they are below the inclusion limits of the catalogue, which are a 8a detection in polarized intensity and P21 > 0.5%. For 
these sources also a significant redshift evolution can be seen with a significance level of 99.64%. 
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Fig. 3. — Comparison of the cumulative RM distributions of BOS 
(solid) and TSS09 (dash) for QSO sight lines with no Mgll ab- 
sorption systems (black lines) or with one or two, respectively (red 
lines). RM data in TSS09 were available for 54 out of the 71 QSOs 
in B08. 



is the theoretical value of the degree of polarization 
which depends on the spectral index 7 of the emit- 
ting relativistic electrons (~ 0.74 for 7 ~ 2.8), and 



Po is the degree of polarization as A — > 0. 

ii) A further effect which reduces the degree 
of polarization is d ifferential Faraday Rotation 
(|Sokoloff et alJlT99l . This effect happ ens if the 
synchrotron emission and the medium producing 
the Faraday rotation are not spatially separated. 
In such a case polarization angles from different 
depths within the source are rotated differently. 
The line of sight integrated emission suffers increas- 
ing depolarization with increasing wavelength. For 
extragal actic sources t his effect seems not to be im- 
portant (|Tribbldll99il) . 

iii) Inhomogeneous Faraday Rotation screens within 
the radio beam lead to depolarization of the 
sources. If the RM screen is modelled by many 
independent RM cells, this effect is cal led depolar- 
ization by external Faraday dispersion (|Burnlll966l : 
ISokoloff et aiiri998l ). If the RM screen varies sys- 
tematically within the radio beam this effect is 
called beam depolarization. 

iBuTLI (fl966l) gives the formula 

p = Pa cxp(-2er|j M A 4 ) (3) 

to describe the depolarization induced by an inhomoge- 
neous Faraday screen with a RM dispersion <t_rm- As- 
suming for simplicity that each RM cell is a cube of 
linear size Iq, then each cell contributes a dispersion 
a c = 0.8lBnJo and ctrm = < J c^/v, where r\ — L/Iq is 
the number of cells along the path-length L traversed by 
the radio waves. 
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The size of the cell is determined by the scale above 
which the RM contributions are uncorrelated and can be 
determined by computing the structure function in high 



resolution RM mapsQ 
IRossetti et all (2008 



(|2008D and IMantovani et all 1)20090 
modelled polarisation observations of ~ 65 compact 
steep spectrum sources between 2.8 to 21 cm done with 
the WSRT, VLA and Effelsberg telescope. They ob- 
served that the fractional degree of polarisation at large 
wavelengths is too large to be explained by Burn's depo- 
larization law (Eq[3]). They observed that for a large frac- 
tion of sources p remains approximately constant above 6 
- 13 cm. To account for these observation they suggested 
that just a fraction of the polarised source is covered by 
a depolarising inhomogeneous RM screen and modified 
Eq.Elto 



P = Po(/ c exp(- 



-24 M A 4 ) + (1 - / c )), 



(4) 



where f c is the covering factor of the source. 

Here we emphasize that partial coverage of the polar- 
ized source by inhomogeneous RM screens is the key for 
explaining the differences in the RM values of K08 and 
TSS09. 

4. MODELLING 

4.1. Extension to cosmological screens 

In Eq. [4] the (1 + z )~ 2 correction for cosmological dis- 
tances is not included. Assuming that the RM dispersion 
causing the depolarization, orm, is constant with z, we 
modify Eq. 2] for cosmological sources as: 

p =p (/cexp(-2a| M (l + z)~ 4 A 4 ) + (1 - f c )). (5) 
Thus the width of p(X) changes as a function of z as 



(1 + ^) 2 

2o-RM 



(6) 



This shows that depolarization by a non-evolving rest 
frame Faraday dispersion screen is expected to decrease 
for screens at higher redshift. Below, we explicitly com- 
pute the depolarization by typical RM screens as a func- 
tion of redshift of the sources. 

4.2. A simple model 

In our model the Faraday screens consist of three com- 
ponents: (i) an inhomogeneous foreground screen local to 
the source with covering factor f c ,QSOi (ii) an inhomo- 
geneous foreground screen in intervening galaxies with 
covering factors f c ,Mgii and (iii) a homogeneous screen 
due to the Milky Way, which is assumed uniform across 
the extension on the sky of the polarized source. 

The uniform Faraday screen in the Milky Way is set 
to a constant val ue, RMmw = —1 rad m 2 . This value 
is consistent with iSchnitzelen (|2010f ) who determined the 
Milky Way c ontribution o f Rota tion Measures at \b\ > 
20° using the lTavlor et all 1)20091) data. 

On the other hand, each inhomogeneous foreground 
Faraday screen is characterized by an ensemble of cells of 
size lo with uncorrelated RM values. The RM values have 

1 The structure function D at scale s is denned as D(s) = 
((RM(x + s) — RM(x)) 2 ) and (} means ensemble averaging. 



a Gaussian distribution with dispersion ctrm^AI + z x ) 2 
and zero mean, where a x is the rest-frame RM dispersion, 
z x the screen's redshift and, x — QSO, Mgll, labels dif- 
ferent screens types, for which the above parameters may 
differ. 

If the RM screen is at the source redshift and the source 
has linear size s, then it is covered by N = s 2 /Iq cells. Of 
these, only a fraction / CjS will be Faraday active for each 
screen. If the RM screen is at the redshift of the Mgll 
systems, s is the projected linear size of the source viewed 
from the Earth at the distance to the Mgll system. 

Each screen is then fully described by a realization of 
N RM values, a fraction f CyX of which are extracted from 
the associated Gaussian distribution to characterize the 
active cells, while the reminder are set to zero to rep- 
resent the inactive cells. Assuming a uniform flux of 
normalized intensity (1=1) and uniform (zero) intrinsic 
polarization angle across the source surface, we can write 
the following relations for the remaining non-zero Stokes 
parameter U, Q, the angle and degree of polarization, 
respectively, 



U(X 2 ) = ^Y,sm(2RM. l X 2 ) 
1 N 

Q(X 2 ) = — ^ cos (2RM t X 2 ) 



X(A 2 ) 
P(A 2 ) 



i=l 

1 fu 

— arctan — 

2 \Q 

VU 2 + Q 2 , 



where 



RMi — RMgso,i + RMM g ii,i + RMmw, 



(8) 

(9) 
(10) 

(11) 

are random variables determined by a Monte-Carlo real- 
ization. 

In Figure 0] we present generic results for three rep- 
resentative cases, namely a low {zqso < 1) and high 
{zqso > 1) redshift QSO with no intervening absorber 
in the top and mid panel, respectively, and a high red- 
shift QSO with intervening absorber in the bottom panel. 
More detailes are specified below. In all cases, the inho- 
mogeneous screen comprise eight cells, N f c>x = 8. 

In each panel we plot the polarization angle, x (black 
solid line - left y-axis), and the polarization fraction, 
p (red solid line - right y-axis) as a function of A 2 . To 
emphasize the importance of the effects of depolarization, 
we also show the rotation in polarization angle produced 
by (a) the average RM within the beam (dash line) , 



X = RM avg X 2 , 



where 



1=1 

1 N 

=RM MW + JjY. RM QSO,i + RM Mg ii,i, (13) 

and (b) by the Milky Way (dash-dot line), 
X = RM MW X*. 



() 



Bernet ct. al. 



Note that RM avg is in general non-zero and different 
from RMmw ■ This is because for fixed fraction of active 
cells the second term in Eq. [13] is a random number of 
aplitude 1/yN, with null probability of being either zero 
or — RMmw- Finally, from left to right, the vertical blue 
dash lines indi cate the typical observational wavelengt h 
of 20 GHz for (IJackson et al.ll20Tot IMurphv et aDl201(fl . 
of 5 GHz for K08 and of 1.435 GHz and 1 365 GHz for 
the NV SS survey (jCondon et al.lll998l )) and lTavlor et~all 
(2003). 

The top panel illustrates the case of a low redshift 
sources (z < 1.0), with observed RM dispersion of 
<trm,QSo/{1 + z) 2 = 28 radm" 2 , f c , Q so = 0.5 and no 
intervening Mgll absorber (f c ,M g ii = 0). It can be seen 
that at short wavelengths, A 2 < 0.014, application of 
the simple A 2 relation in Eq. [TJ yields the average RM 
within the beam from which the R M distributions in th e 
foreground screens can be inferred (jBernet et al.ll2008h . 
However, at longer wavelengths the flux through the in- 
homogeneous screen is depolarized so that the change in 
polarization angle is dominated by the contribution of 
the Milky Way. Correspondingly, the degree of polariza- 
tion is given by p ~ 1 - f c ,QSO- 

The mid panel illustrates the case for a high red- 
shift source, z > 1.0 with observed RM dispersion of 
&ralqso/(1 + z) 2 = 8 radm" 2 and J c ,qso = 0.5. In 
this case, depolarization is negligible up to A 2 = 0.05 or 
at 1.4 GHz. This means across most of the considered 
wavelength range application of the A 2 relation in Eq. 
Q] would yield the average RM contributed by both the 
QSO screen and the Milky Way. 

Finally, the bottom panel illustrates the same case as 
the mid panel but with an additional Faraday screen due 
to an absorber, with <jRM,Mgii /(l + z) 2 ~ 230 rad m -2 
and fc,M g ii = 0.5. 

This case has direct connection to the study of B08, 
who showed that intervening galaxies traced by Mgll ab- 
sorption lines contribute an additional Faraday screen 
with an observed rest frame RM dispersion OMgii = 
140 radm~ 2 . In the context of our simple model, this 
implies a RM dispersion of the individual cells of the 
screen, responsible for the depolarization, <JRM,Mgii = 
CMgiiy/N/ fcAigii- The value for URMMgii/ (1 + z) 2 as- 
sumed above is is based on this relation and the choice 
of N f c ,Mgii = 8. The bottom panel of Fig. |4] shows 
that with this choice of parameters the intervening screen 
causes significant depolarization above A 2 ~ 0.003m 2 . In 
this case the application of the A 2 relation at longer wave- 
lengths yields the Milky Way RM value. 

The approach taken here, is valid in the limit of large 
number of random cells. U sing the total flux - an gular 
diameter relation given by iWindhorst et all (|1984j i we 
can estimate the extension of the region of the polarized 
emi ssion. The median flux of our sources which are in 
the lTavlor et al.l (|2009( ) catalogue (371 sources) is 1.2 Jy. 
This translates to a median source size in total flux of ~ 
10 arcsecs. Since the polarized emission is likely to be 
dominated by some hotspots the effective size in polar- 
ized emission is expected to be smaller than that. Given 
that the extension of the polarized emission is a few arc- 
secs and what we know about the scale of fluctuations of 
the RM (see next section) we expect a significant number 
of cells within the beam. In addition, using Eq. [5] and 



future multi-wavelength radio polarization observations 
one can determine cr RM,Mgii and f c ,Mgii directly and 
thus calculate N. 

Finally, we hav e assumed a homogene ous RM screen in 
the Milky Way. IGaensler et al.l ([2005D . however, shows 
that arcsec scale fluctuation exi st in the RM sc r een o f 
the Large Magellanic Cloud and Havcrko rn et al.l (2008) 
shows that there are RM fluctuation in the local inter- 
stellar medium. Inhomogeneous Galactic RM contribu- 
tion, which might vary for different sight lines, would 
lead to additional depolarization and add scatter to the 
predicted trends (see section 14. 3|) of our simple model, 
but cannot mimic them. 

4.3. Model predictions 

Based on this toy model with the free parameters f c ,x, 
o~rm,x, we can predict simple trends in the observed angle 
and degree of polarization measured at fixed frequency 
as a function of radio source's redshift. For simplicity we 
shall assume there is no redshift evolution in the model 
parameters. We must also distinguish between the two 
cases with and without the Faraday screen provided by 
intervening galaxies. 

For the case without intervening galaxies the model 
predicts 

i) that towards low redshifts P2i/p\^a ~ 1 — f c ,QSO, 
because the flux through the inhomogeneous Fara- 
day screen is depolarized. 

ii) that towards high redshift P2i/p\^o ~ 1 because 
the depolarization effects are suppressed. How- 
ever, P2i/pa->o anti-correlates with ctrm,qso due 
to residual depolarization effects associated with 
the local screen. 

iii) that towards low redshifts the discrepancy 
|RMa^o — RM21 1 increases because due to depolar- 
ization effects RMa-j-o measures RM avg while RM21 
is dominated by the Milky Way contribution. 

iv) that the discrepancy |RMa^o — RM21I decreases 
towards high redshifts because both RM^_> and 
RM 2 i tend to measure the average RM avg (Eq.lT2"j). 

For the case with intervening galaxies the model predicts 

i) that at all redshift the discrepancy 
|RMa^o — RM 2 i| is large, because due to de- 
polarization effects RMa->o measures RM avg and 
RM 2 i is dominated by the Milky Way contribution. 

ii) that at all redshift P2i/pa^o 1 — fc,Mgii because 
the sources will be significantly depolarized. 



5. DISCUSSION 

5.1. Nature of the partial inhomogeneous RM screen 

It is instructive to look at high resolution Very Long 
Baseline Array (VLBA) radio polarization observations 
of some of our objects in K08. Examples of sources 
where VLBA polar ization o bservations exist are 3C43 
(ICotton et al.H2003D. 3C118 (IMantovani et al.l l200ll and 
B1524-026 (jMantovani et al.ll2002D . Typical resolutions 
of these observations are ~ 8 x 8 mas 2 , which corresponds 
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Fig. 4. — Angle of polarization, \ (black solid line - left y- 
axis), and theoretical polarization fraction, p (red solid line - right 
y-axis), as a function of A 2 , for three models of inhomogcncous 
Faraday screens: a low (top) and high (mid) redshift QSO with 
no intervening absorbing system, and a high redshift QSO with 
intervening absorber (bottom). The dash line indicates the rotation 
of polarization angles due to the average RM within the beam, x = 
RM avg X 2 , and the dash-dot line that due to the Milky Way, \ = 
RMmw^ 2 ■ The vertical blue dash lines from left to right, indicate 
typical wavelength of observations at 1.5 c m, at 6 cm (K08), an d 
at 20.89 cm and 21.96 cm (NVSS survey, ICondon et all (Tl998h ) 
for lTavlor et alJ l|2009l) calculated RM values. In the top and mid 
panel the Faraday screens include an inhomogeneous screen local 
to the source with a covering factor f c = 0.5 and a screen in the 
Milky Way that is uniform across the source extension. In the top 
panel the source QSOs is at z < 1.0 and has ctrm,QSo/(1 + z) 2 = 
28 radm" 2 ; in the mid panel the QSO is at z > 1.0 and has 
&RM,QSO /(I + z ) 2 = 8 radm -2 . The bottom panel is the same 
as the mid panel with an additional Faraday screen due to an 
intervening system, with ctrm Mgll /(I + z ) 2 ~ 230 radm" 2 and 
/ c = 0.5. 



to a resolution of 64 pc x 64 pc at z ~ 1. At such a 
high resolution complex structures in polarization angles 
and RM m aps are revealed. Ofte n it can be seen, e.g 
B 1524-026 ([Mantovani et al.ll2002h that there is a dom- 
inant compact component in polarised flux and a more 
extended diffuse polarised component. Further it can 
be seen that the diffuse component consist of many in- 
dependent RM cells. That means that for unresolved 
observations at large wavelengths the diffuse component 
will cancel and the more compact component dominates 
the observations. 

For some sources it possible that the number of cells 
f c N in the RM screen is very low. In this case no depolar- 
ization is observed but p and x oscill ate. This situation 
was observed for the source 3C 27 bv lGoldstein fc Reedl 
(|1984fl . iRossetti etail (|2008| ) fitted a two component 
model to the data to describe p and \ for the source 
B3 0110+401. iFarnsworth etaLl (poTl also used a two 
component model to describe radio polarization observa- 
tion at large wavelengths A~lm. 

For the sightlines with intervening galaxies it is very 
natural to assume that the magnetic fields within them 
wil l lead to depolarizatio n, an effect studied in the LMC 
bv lGaensler et all (|2005D . Using the Milky Way as typ- 
ical galaxy we would expect depolarization of the back- 
gro und QSOs. For NGC 1310 t his effect was observed 
by iFomalont et aTl (|1989t ) and iSchulman fc Fomalontl 
(|1992l ). One possibility is that the source covers both the 
spiral arms and interarm regions. For the spiral arms the 
coherence lengths of the magnet ic fields are much shorte r 
as suggested by observations of lHaverkorn et ail (|2008|) . 
These authors gives an outer scale of turbulence of 0.1 
kpc for interarm regions and 10 pc for the spiral arms 
in the Milky Wajo In this way the area covered by 
the spiral arms would be depolarized and the interarms 
would produce a coherent screen. Other interpretations 
are again as for the QSO itself that the polarized emis- 
sion is dominated by one compact component which com- 
petes with a more extended diffuse emission. If the size 
of the compact component is small enough, this compo- 
nent could get a coherent screen whereas the extended 
source gets different RM screens which will lead to depo- 
larization of this component. 

5.2. Rotation Measure Synthesis 

At the moment there are several pol arization surveys 
under way or plan ned, e.g. GALFACTS dTavlor fc Salterl 
[2010I). LOFAR (lAnderson et al.1 [2011 . or POSSUM 
(|Gaensler et all l2010f ) within ASKAP. These polariza- 
tion surveys should pave the way for the planned polar- 
ization survey for SKA dGaensler et al.1 [20041 : iGaenslerl 
2009J iBeck fc Gaenslerl 12004ft . In particular one wishes 
to study the evolution of magnetic fields in galaxies and 
in the intergalactic medium as a function of cosmic time 
using large RM datasets. 

With the large number of available spectral channels 
in these surveys, e.g. a few thousands, Faraday Ro- 
tation Measure Synthesis (jBrentiens fc de BruvrJ 12005ft 
can be performed. Using this technique one performs a 

2 The small coherence lengths of the magnetic fields has probably 
to do with our point of view in the disc of the Milky Way. We know 
of course from observations of nearby spiral galaxies that there is 
a magnetic field component which is coherent over several kpc. 
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Fig. 5. — Original Faraday spectrum for a low redshift QSO 
and its associated observed polarization angle and degree, together 
with the recovered Faraday spectrum for a finite wavelength cov- 
erage. The Monte Carlo realization is identical to the top panel 
in Figure [4] Upper left panel: Original Faraday spectrum for a 
low redshift QSO with a homogeneous Milky Way RM screen of 
RM = — lOrad m 2 and an inhomogencous intrinsic RM screen with 
a RM,QSo/(l + z ) 2 = 28 rad m -2 and a covering factor f c = 0.5. 
Upper middle panel: Corresponding U and Q Stokes parame- 
ters (thin solid line, thin dashed line) and degree of polarization 
p (thick solid line) as a function of A 2 . An observational error 
a U,Q ^ 0.017% is added to every channel. Upper right panel: Cor- 
responding polarization angle x vs. A 2 . The dashed line indicates 
the expected polarization angle from the homogeneous Milky Way 
RM contribution. Lower left panel: Recovered Faraday spectrum 
for the wavelength range X m i n = 0.03m to X ma x = 0.27m. The 
thin solid/dashed line indicates the real/imaginary part of the spec- 
trum and the thick solid line shows the amplitude. Lower middle 
panel: Faraday spectrum after applying the CS-RM-Thin method 



panel: faraaay spectrum alter applying tne os-ruvl-ll 
of lLi et al.l 1120111 ) to correct for the limited wavelength 
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Fig. 6. — Identical to Figure \E\ but for the case of a high redshift 
QSO with an intervening Mgll system. The Monte Carlo realiza- 
tion is identical to the bottom panel in Figure [4] 



Fourier transformation of the observed complex polar- 
ization p(A 2 ) to obtain a Faraday depth spectrum F(4>). 
Here the case of a single uniform Faraday screen corre- 
sponds to a delta function in Faraday depth space and 
in this case the Faraday depth is equal to the traditional 
RM(Eq.p. 

Basic parameters for a Faraday survey are the max- 
imum observable Faraday depth, 4> max ~ V3/S\ 2 : the 
Faraday depth resolution, <5</> w 2-^/3/ AA 2 and the largest 
scale in Faraday depth that one is sensitive to, max-scale 
7r/A^j in . Here <5A 2 is the channel width, A A 2 is the 
total bandwidth and X min is the smallest wavelength 
(squared) of the observations. Regions which produce 
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Fig. 7. — Identical to Figure[6]with same AA 2 coverage and thus 
identical Faraday depth resolution Stj) but with observations done 
at longer wavelengths. 



only Faraday rotation but do not emit polarized radiation 
can be described by D irac S functions in (j) space an d are 
Faraday-thin sources (jBrentiens fc de Bruvni r2005'). Be- 
low we discuss our model of the depolarization of QSOs 
in terms of Faraday Rotation Measure Synthesis. This is 
to provide some basi c comparisons for ongoi ng or future 
RM surveys. See e.g. lO'Sullivan et al.l (|2012l ) for current 
observations of Faraday depth spectra. 

In the next, with the auxilium of Fig. <(5])— ([T]) we give 
a simple example of the application of RM-synthesis to 
the same cases studied in Figure 21 which illustrates the 
advantages of this emerging technique. In Figure [5] the 
original Faraday spectrum for a low redshift QSO with- 
out intervening galaxies is shown. Here the Monte Carlo 
realization is the same as in Figure |4] (upper panel) with 
adopted parameters urm,qso/{^ + z) = 28 radm -2 
and fc,QSO = 0.5. In particular, all of the 16 cells have 
the contribution from the MW and eight cells have an 
additional intrinsic low-z QSO contribution. The Fara- 
day spectrum is only real valued which means that the 
intrinsic polarization angles are zero. The adopted mini- 
mum and maximum wavelengths are Xmin = 0.03 m and 
Xmax = 0.27 m and correspond to a Faraday depth res- 
olution of Sep « 48.1 radm -2 . The choosen wavelength 
range corresponds to the offered wavelength range at the 
Australia T elescope Compact Array which is similar to 
the EVLA (|Beck et al.l 12012ft . The recovered Faraday 
spectrum with this resolution is shown in the lower left 
panel. 

There are several methods to recover the informa- 
tion lost due to the incomplet e wavelength coverage 
IjHeald et al.l(200llLi et alJl20Tl . We use here the com- 
pressive sampling (CS) method of iLi et all (|2011ft . Best 
results are obtained if prior information about the Fara- 
day spectrum is present, e.g. if the sources are Faraday 
thin or thick. The obtained Fara day spect r um af ter ap- 
plying the CS RM thin method of ILi et al.l (|2011| ) to the 
Q and U vectors (upper middle panel) is shown in the 
lower right panel. It can be seen that the homogeneous 
RM component at -10 rad m -2 can be approximately re- 
covered, but not all of the inhomogeneous RM compo- 
nents. 

In Figure [S] the Faraday spectrum for a high redshift 
QSO with an intervening galaxy is shown. Here the Fara- 
day spectrum is identical to the one used in Figure [5] in 
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the lower panel with the parameters <JRM,Mgii / (^+z) 2 ~ 
230 radm -2 and f c<M gll = 0.5, a RM ,QSo/{^ + z) 2 ~ 
8 radm -2 and f c ,Qso — 0.5, and RM = — 10 radm -2 . 
In the Monte Carlo realization shown in Figure |5] all of 
the 16 cells have the contribution from the MW, eight 
(random) cells have an additional intrinsic high-z QSO 
contribution and eight (random) cells have an additional 
contribution from the intervening galaxy. With the ob- 
servational parameters \ m in = 0.03m and X ma x = 0.27m 
and 5(j) ps 48.1 radm -2 the MW+QSO components can 
not be resolved but lead to a distinct peak at <f> rs 
— 10 rad m -2 with an amplitude ~ 0.45 m rad -1 in 
the recovered Faraday spectrum. On the other hand, 
the eight RM components with contributions from inter- 
vening galaxies can be well resolved with the assumed 
covered wavelength range. 

For comparison in Figure [7] we show the recovered 
Faraday spectrum for observations using longer wave- 
lengths, Xrnin = 0.27 m and X max — 0.38 m but with the 
same covered wavelength range AA 2 . This shows that 
the qualitative features of the Faraday spectrum are also 
recovered, although with a lower quality than in Fig. [6] 

In summary, from section 14.31 using polarization angle 
X vs. A 2 observations one might conclude that it is best 
to do observations at short wavelengths (below the ex- 
ponential fall off) in order to measure the average RM. 
However as we have illustrated above using RM Synthe- 
sis this is not necessary. In order to be able to resolve the 
individual RM co mponents a large AA 2 is required. (See 
IBeck et al.l IpOll for a summary of the Faraday depth 
resolution 8<j) of current and future radio telescopes.) The 
proposed wavelength range for SKA with A ~ 0.03— 4.3m 
and AA ~ 18m 2 and 5<f> ~ 0.2 radm -2 will offer a superb 
Faraday resolution to resolve individual RM components 
in intervening galaxies. 

6. SUMMARY AND CONCLUSIONS 

In previous work we have used RM values from K08, 
typically determined at 6 cm, to study magnetic fields in 
normal galaxies at high redshift. When the analysis was 
repeated using the recently released RM data of TSS09, 
determined at 21 cm, we could find neither a correlation 
between |RM| and absorbing systems as in B08, nor an 
increase in the RM dispersion with z as found in K08. 
Motivated by the above results, we have attempted to 
understand those differences in terms of a simple model 
based on inhomogeneous Faraday screens associated both 
with the QSO and the Mgll host galaxies, and for both 
low and high redshift. 

We find that the presence of inhomogeneous screens 
leads to important departures from the classical A 2 de- 
pendence of the rotation of the polarization angle. Re- 
lated to this are depolarization effects which become 
stronger towards higher wavelengths. As a result, due to 
cosmological redshift, observations at fixed wavelength 
are affected differently by depolarization and are sensi- 
tive to different Faraday active components. In particu- 
lar depolarization effects become stronger for lower red- 
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shift QSOs. We find that the depolarization saturates 
to a value given by the intrinsic polarization times the 
complement of the covering factor of the inhomogeneous 
screen, (1 — f c ). 

Actual predictions depend on the assumed values for 
the model parameters. The following results apply for 
the choices made in section. 21 which are relevant for 
the current investigation. When the line of sight to the 
QSO is free of absorption systems, application of the A 2 
analysis to extract RM values from radio observations 
then has the following consequences: 

• for low redshift QSOs, RM21 is dominated by the 
Milky Way contribution, while RM^o measures 
RM avg (see Eq. H5J. 

• for high redshift QSOs, in general the RM value 
reflect RM avg 

The presence of absorption systems contributes an ad- 
ditional Faraday screen which further depolarizes the low 
frequency radiation. Therefore, 

• while RMa^o is given by RM avg , RM 2 i is domi- 
nated by the Milky Way contribution and the dis- 
crepancy |RMa^o — RM21I is large at all redshifts. 

In conclusion, while the model is admittedly simple, it 
seems plausible to consider that the discrepancy between 
results based on K08 and TSS08 RM are due to the se- 
vere depolarization induced by inhomogeneous Faraday 
screen on high wavelength radiation. This, however, does 
not exclude the importance of other effects. 

Finally, the study of magnetic fields in galaxies as a 
function of cosmic time will benefit considerably from the 
application of RM-synthesis, which has the power to dis- 
intangle the contribution from inhomogeneous magneto- 
active components. For this purpose, instrumental band- 
width is most important, although higher frequency 
channels appear to deliver higher quality. This conclu- 
sions may be refined with future investigations. 

We are very grateful to an anonymous referee for 
several valuable comments which helped improve the 
manuscript considerably. This work was supported by 
the Swiss National Science Foundation and has made use 
of NASA's Astrophysics Data System. During the refer - 
eeing process a related work bvlArshakia n fc Beck! f|201 lh 
appeared. It makes general predictions for future radio 
polarization surveys in qualitative agreement with ours, 
though it does not study the case of intervening gala xies. 
Further the RM catalogue of lHammond et ahl (|2012t ) ap- 
peared with determ ined redshifts of 4003 sources of the 
iTavlor et ail (|2009f) sample. Similar to our findings they 
do not see any redshift evolution of the RM distribution 
in their large sample. 
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